Workshop on the CKM Unitarity Triangle, IPPP Durham, April 2003 



o 
o 

(N 
C 

(N 



> 

O 

m 
o 

a: 

Oh! 
<D 
43 



X 



Prospects for B d s —> h + h measurements at LHCb 

V. M. Vagnoni 

Dipartimento di Fisica di Bologna and INFN, Sezione di Bologna, via Irnerio 46, 1-40126 Bologna, Italy 

LHCb will collect large samples of B d and B s decays. Combining the CP-violating observables of the decays B d — > n + n~ and B s — > K + K~ 
it is possible to extract the y angle of the unitarity triangle. The selection of these decays within the current LHCb simulation framework 
is outlined and the expected annual event yields and background-to-signal ratios are quoted. Then, the results of a study on the sensitivity 
that LHCb can achieve for the corresponding CP-violating observables are presented. 




Figure 1. Schematic drawing of the LHCb detector. 



1 Introduction 

LHCb is a dedicated experiment on b-quark physics, cur- 
rently under construction at the Large Hadron Collider 
(LHC). It will profit from the 500 pb b-hadron produc- 
tion cross section in the 14 TeV proton-proton collisions 
at LHC (i.e. about 1% of the total visible cross section) 
to make precise measurements of CP-violation and rare 
decays of the B-mesons [QJ. By over-constraining the 
Cabibbo-Kobayashi-Maskawa matrix elements, LHCb will 
hopefully be able to observe subtle inconsistencies with 
the Standard Model, therefore providing indications of new 
physics. A sketch of the LHCb detector is shown in Fig. ^ 

Profiting from the large bb cross section, LHCb will collect 
large samples of Bd and B s decays. Combining the CP- 
violating observables of the decays Bd — > tt + tt~ and B s — > 
K + K~ it is possible, by means of a method proposed by R. 
Fleischer [ 2 1, to extract the y angle of the unitarity triangle. 

In section 2 the selection strategy, implemented to iden- 
tify the Bd — > n + n~ and B s — > K + K~ decays within the 
current LHCb Monte Carlo simulation framework, is de- 
scribed and the estimated annual event yields and back- 
ground levels are quoted. 

By using the results shown in section 2, the sensitivities 



Channel 


Monte Carlo Statistics 


B d -» 7T+7T 


60000 


B s -> K + K~ 


65000 


B d -> K + n- 


62000 


B s -» K n + 


24000 


A h -> pK~ 


20000 


A b -> piC 


21000 


bb-^X 


1144000 



Table 1. Monte Carlo event samples used to study the selection 
of B,j — > n + n~ and B s — > K + K~ decays. The other channels are 
studied as potential sources of background. 



achievable on the CP-violating observables can be deter- 
mined. This is described in section 3. 



2 Event reconstruction 

To evaluate the performance of LHCb in the reconstruction 
of B d -> n + n- and B s -> K + K decays, a full GEANT 
Monte Carlo simulation has been performed, followed by 
realistic pattern recognition algorithms for the track re- 
construction and particle identification. The Monte Carlo 
event samples used for this analysis are shown in Tab. ^ 

Two sources of background have been considered: two 
charged body decays of B mesons and baryons, which 
can fake the signal in case of particle mis-identification, 
and combinatorial background. The dominant source of 
combinatorial background is believed to come from beauty 
events. 

To select the interesting decays and reject the backgrounds, 
a set of kinematical and topological cuts has been used. For 
each pair of tracks with opposite charge identified as pions 
or kaons by the RICH detectors, cuts are applied on: 

• momentum; 

• smallest and largest transverse momentum; 
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Figure 2. Invariant mass distribution for B d — » n + n~ recon- 
structed events, together with the mass spectrum due to other 
two body B-decays background. The mass resolution is about 
ISMeV/c 2 . 



• smallest and largest impact parameter significance; 

• x 1 of common vertex fit. 

Each pair surviving to these cuts is used to form a B-meson 
candidate, and further cuts are applied on it: 

• transverse momentum; 

• impact parameter significance; 

• distance of flight significance; 

• invariant mass. 

With these set of cuts, with values specifically optimized 
for each of the two decays under study, all the combina- 
torial bb background events, in the limited Monte Carlo 
statistics available at the moment of this analysis, are re- 
jected. At the same time the two body B decays studied as 
specific backgrounds, thanks to the particle identification 
performance relying on the two RICH detectors [ |3, are 
maintained under control. 

Fig. 13 shows the invariant mass distribution for the re- 
constructed Bd — > 7r + 7r~ events surviving the selection, to- 
gether with the two body B-decay background. The esti- 
mated annual event yields and background-to-signal ratios 
can be read in Tab. [2] 

3 CP sensitivity 

The time dependent CP asymmetry takes the usual form: 



9\cp(t) = 3K% ■ cos(x ■ t) + $&Zsin{x ■ t) 



Channel 


B/S 


Untagged annual yield 


B d -» 7T+7T 


< 0.8 


27000 


B s -> K + R- 


<0.55 


35000 



Table 2. Estimated number of offline-selected events per year and 
background-to-signal ratios for the B d — > n + n~~ and B s — > K + K~ 
decays at LHCb. 



which is valid assuming the decay width difference of the 
B mass eigenstates AT to be negligible. 

In order to generate the requested samples of tagged Bd —> 
jt + jt~ and B s — » K + K~ decays to perform a study on the 
sensitivity that LHCb can achieve on the measurements 
of and JV^'p, a standalone toy Monte Carlo program 
has been used. The program takes as input the number 
of reconstructed events per year and the background levels 
shown in Tab. |3 It reproduces the acceptance as a func- 
tion of the proper time and the proper time resolution as 
studied in detail by using the full GEANT simulation and 
simulates the effect of CP violation using given values of 
Jvqp and M™p, then simulates the tagging procedure, fi- 
nally giving as output a sample of proper times of tagged 
B and B decays. 



To "measure" from this sample of proper times the values 
of J?lip and J?Kp , with their uncertainties and correlation, 



an unbinned maximum likelihood method has been used. 
The log-likelihood function to be maximised with respect 
is: 



to J[f p and W?' x 



'CP 



log Z = log 



[~[h(t,)-]~[«(t,) 



where {r,} and {r,} are the sets of decay proper times for 
the tagged B and B respectively, and 



n(T;) = 



n(Tj) 



,-1/tb 



\-2oj \ 

3\cp(t)\- e(t) 

1 + BIS 



Rto - 1) 



R(Tj - t) 



being Tp the B-meson lifetime, B/S the background-to- 
signal ratio, to the wrong tagging fraction, e(t) the accep- 
tance as function of the proper time, R(t - 1) a function ac- 
counting for the proper time resolution of the LHCb spec- 
trometer (normal distribution with a width of 40 fs) and 
where the symbol stands for convolution product. 

As an example, Fig. [5] shows the dependence of log £, on 
3Tq P and for a sample of Bd — * n^iC and a sam- 
ple of B s — > K + K~ events generated by the toy Monte 
Carlo program, corresponding to one year of LHCb data 
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Figure 3. Dependence of the log-likelihood on and !K™p for 
the Bj — > n^n" decay (upper plot) and for the B s — » K + K~ decay 
(lower plot). 



Channel 








B d -> tt+tT 


0.054 


0.054 


-0.53 


B s -> K + K- 


0.043 


0.043 






Table 3. Resolutions on the measurements of M^' p and tftgp for 
the B c \ — > n*n~ and B s — > K + K~ decays expected at LHCb after 
one year of data taking. The last column reports the correlation 
between the measured values of the asymmetries. 

taking, where the following input values for the CP asym- 
metries have been used: JV^'piB^ — » n + n~) = -0.3, 
y^{B d -> tt+tT) = 0.58, mf p {B s -> K + K ) = 0.16 
and 3K™(B S -> K + K ) = -0.17. Fig. E| shows the cor- 
responding confidence regions determined for y\^ r p and 
y[™p , while Fig. [5] shows the binned asymmetry for the 
Bd — » 7r + 7r~ decay with the result of the likelihood fit su- 
perimposed. 

The obtained resolutions and correlations for &Hj?p and 
y^Qp, corresponding to one year of LHCb data taking, are 
shown in Tab. [3] The measurements of 3y£p and 9Vgp are 
strongly correlated in the case of the B d — > n + n~ channel, 
while the correlation for the B s — » K + K~ channel is negli- 
gible, as it can easily be argued by looking at the contour 
plots of Fig. |U 
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Figure 4. Confidence regions for and for the B d — > 
7r + 7r~ decay (upper plot) and for the B s — > K + K~ decay (lower 
plot), corresponding to one year of LHCb data taking. The four 
ellipses correspond to the 39% (standard error ellipse), 68%, 95% 
and 99% confidence regions. 
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Figure 5. Binned asymmetry for the Bj — > n + n decay with the 
result of the likelihood fit superimposed. 
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Figure 6. Dependence of the error on && r p for the B s 
decay on the value of the mixing parameter x s , 



K + K- 



The previous results for the B s — > K + K~ channel are ob- 
tained for a value of the mixing parameter x s = 20. Due to 
the finite proper time resolution of the spectrometer (about 
40 fs), higher values of x s , i.e. faster B s - B s oscilla- 
tions, lead to larger errors for the measurements of 
and tt™p . Thus, a study of the dependence of the reso- 
lution of SflPgp and JVg* on x s , up to x s = 40, has been 
performed. The results of this study for Jl^p are shown in 
Fig. [6] (the dependence for JV™p is identical and is omit- 
ted here). As it can be seen the error on the asymmetries 
increases by about a factor 1 .6 at x s = 40 with respect to 
x s = 20. 



4 Conclusions 

The results presented in this paper show the potential of 
LHCb in collecting large samples of Bj — » 7r + 7r~ and 
B s — > K + K decays. They are obtained by using a full and 
realistic Monte Carlo simulation of the detector response. 

The number of Bj — > n + n~ offline-selected decays are 
27000/year, while for B s K + K it is 35000/year, with 
background-to-signal ratios respectively of about 0.8 and 
0.55. 

The sensitivity on the CP-violating observables achievable 
by LHCb for these channels has been studied in detail, 
yielding very encouraging results. The estimated precision 
on Ji^p and ZRFgp after one year of data taking, assuming 
x s to be equal to 20, is 0.054 for the B c / — > n + n~ channel 
and 0.043 for the B s -> K + K~ channel. 



By means of these observables it is possible to extract the 
y angle of the unitarity triangle. In fact, J^p and SV^p, 
related to the Bj — » ji + jt~ decay, can be expressed in the 



framework of the Standard Model as: 
2d sin 6 sin y 



rnimx/ p \ _ 



1 - 2d cos 8 cos y + d 2 

sin(0rf + 2y) - 2d cos 8 sin(<7)</ + y) + d 2 sin <f>^ 
1 - 2d cos 8 cosy + d 2 



where 4>d — 2/3 is the B c i - Bj mixing phase, while d and 
8 parametrize - sloppily speaking - the penguin over tree 
amplitude ratio of the decay transition [|2]. Analogously, 
3K% and related to the B s -> K + K~ decay, can be 

written as 



2d' sin 8' sin y 



1 + 2d' cos 8' cos y + d' 2 



^mix, D v sin (<^ + 2 T) + 2d' cos 8' sm(<t> s 

•SXcpyvs) ~ = 

1 + 2d' cos 8' cos y 



8' sm(4> s +y) + 
+ d' 2 



d' 2 sin <; 



where cf> s = 25y is the B s — B s mixing phase, d' = d' , 

and the d' and 8' parameters are analogous to d and 6 for 
the B s —> K + K transition. 

In the limit of exact U-spin symmetry of the strong inter- 
actions the relations d-d' and 6-8' hold, and the mea- 
surements of the four asymmetry coefficients allow to de- 
termine (f>d and y simultaneously, provided that <p s is deter- 
mined elsewhere (e.g. through the B s — > J/if/<p transition) 
or considered negligibly small, as it is expected in the Stan- 
dard Model. Moreover cpd will be accurately known at the 
time of LHCb, thus allowing a more precise determination 
of y. 

The resolution on y achievable at LHCb by using this 
method, in case of negligible U-spin breaking effects, is 
of the order of few degrees. 
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